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Abstract 

Generalized Dyck paths (or discrete excursions) are one-dimensional 
paths that take their steps in a given finite set S, start and end at height 0, 
and remain at a non-negative height. Bousquet-Melou showed that the 
generating function Ek of excursions of height at most k is of the form 
Fk/Fk+i, where the Fk are polynomials satisfying a linear recurrence re- 
lation. We give a combinatorial interpretation of the polynomials Fk and 
of their recurrence relation using a transfer matrix method. We then ex- 
tend our method to enumerate discrete meanders (or paths that start at 
and remain at a non- negative height, but may end anywhere). Finally, 
we study the particular case where the set S is symmetric and show that 
several simplifications occur. 



1 Introduction and notations 

A Dyck path is a one-dimensional path taking its steps in { — 1,1}, starting and 
ending at and visiting only non- negative points (Figure [TJ left). It is well- 
known that the number of Dyck paths of length 2ri is the nth Catalan number 
C n QUI Chapter 6]: 

_ 1 f2n^ 

Moreover, let E = E(t) be the generating function of Dyck paths. This gener- 
ating function satisfies the following algebraic equation: 



1 — E + t E 



2p2 



0. 



A generalized Dyck path (or discrete excursion) takes its steps in a given 
finite set S instead of { — 1,1}. A discrete meander is a slightly more general 
path: it takes its steps in 5, starts at and visits only non-negative points, but 
may end anywhere (figure [TJ right). 





Figure 1: Left: a Dyck path of height 3. Right: a discrete meander with set of 
steps S = {±1, ±2}, of height 4 and with final height 3. 
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A large number of papers enumerate generalized Dyck paths with varying 
amounts of generality. Methods used include some properties of Laurent series 
[5] , grammars pj [HI [T] , or the kernel method [H [3J . Assume that all steps s 
in S have a weight u s taken in some field of characteristic (typically a field of 
fractions with one or several variables over Q) . Let E and M be the generating 
functions of excursions and meanders, respectively, according to the weights u> s . 
Banderier and Flajolet [2] showed that both these generating functions are al- 
gebraic. More precisely, let a = max S and b = — min S; one may compute 
polynomials of degree ( a ^ b ) canceling E and M. 

In this paper, we consider excursions and meanders with bounded height, that 
is, that never go above a certain level, say k. We denote by E\. the generating 
function of excursions of height at most k and by Ek t i the generating function 
of meanders of height at most k with final height i. We also denote by the 
generating function of meanders of height at most k regardless of final height. 

Bousquet-Melou [5] proved that the generating function Ek is of the form: 

_ F k 
&k — 



+i 



where the Fk are polynomials in the weights u s . She also proved, with the 
use of symmetric functions, that the polynomials Fk satisfy a linear recurrence 
relation of order ( a ~^ b ) ■ In other words, we have 

^ tkZ ~ D{zY 

k>0 y ' 

where the polynomial D(z) has degree ( a ^ b ) (the polynomial N(z) has degree 
(°|') — a — b). Moreover, it can be seen that the polynomial D(z) cancels the 
generating function of excursions E. 

To our knowledge, the only cases of the generating function Ek t t that were 
studied before with general steps are I = 0, that corresponds to excursions, and 
I = k, that corresponds to culminating paths In this paper, we use a transfer 
matrix method to compute both generating functions Ek and Ek,t- 

We now set some notations. As argued in [3], excursions of height at most 
k are walks in a finite graph, with vertices {0, . . . , k} and an arc from i to j if 
j — z is in S. We denote by Ak the adjacency matrix of this graph. If s is in Z, 
let f} s be the quantity: 



ui s if s G S, 
otherwise. 



In this way, the entry of the matrix 1 — Ak is fij-i- 

If m and n are integers, we use the notation [m, n] to denote the set of 
integers i such that m < i < n. We also use the notation N>„ to denote the set 
of integers greater than or equal to n. Finally, if X is a set and n an integer, 
we call n-subset of X a subset of X of cardinality n. 

The paper is organized as follows. Section [2] deals with bounded excursions, 
re-proving Bousquct-Mclou's results using a simple transfer matrix method. 
This method is expanded in Section [3J to cover bounded meanders as well. In 
Section we discuss the case where the set S is symmetric and show that 
several simplifications occur. 



2 



2 Bounded excursions 



Our first step to enumerate excursions of height at most k is the same as in 
[3]. As these excursions are walks that go from to in the graph described 
by the matrix A^, the generating function is the (0,0) entry of the matrix 
(1 — Ak)^ 1 (see [5J Chapter 4]). Therefore, we have: 

Ek = -fi-, (1) 

where Fk is the determinant of the matrix 1 — Ak-i, with the convention Fq = 1. 

As the entry of the matrix 1 — Ak-\ is Pj-i, we have, by the definition 
of the determinant: 

fc-i 

F k = eWlI^H 

rre6 k 1=0 

In the following, we use this expression to compute Fk- The permutations 
a that it involves are, of course, bijections from [0, k — 1] to itself; however, we 
find it more convenient to regard them as bijections from N to itself that fix all 
points in N>k- 

Definition 1. Let I be an a-subset of the set [—6, a — 1]. We call I -permutation 
of order k a bijection from N to I U N> a that fixes all points in N>k- 

Note that for a /-permutation to exist when k < a, all points in [fe, a — 1] 
must be in the set I since they are fixed points. Also note that a [0, a — 1]- 
permutation of order k is the same as a standard permutation of order fc; for 
this reason, we set Jq = [0, a — 1]. 

We denote by ©j^ the set of /-permutations of order k. Let a be a I- 
permutation. We define the number of inversions of a, denoted by inv(er), in 
the same manner as a regular permutation: 

inv(o-) = \ i<j and a(i) > a(j)}. 

We define the signature of a, denoted by e(a), to be the number (— l) mv ( cr ). We 
also define the quantity /3(a) to be: 

fc-i 

/3(a) = e(a) JJ /?<,(,)_,;. 

Finally, we call head of cr the value <r(0); we call tail of a the mapping t defined 
for all n £ N by: 

T (n) = cr(n + 1) - 1. 

We now denote by the vector indexed by the a-subsets of [—6, a — 1] and 
the entries of which are: 

The above remark means that the entry Iq of this vector coincides with Fk ■ 

We compute the vector using a simple transfer matrix method, in a 
manner similar to Example 4.7.7 and Proposition 4.7.8a]. If I is a subset of 
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[—6, a — 1] and s an integer, we denote by e s (I) the number —1 to the power of 
the number of elements of / lower than s: 



ej{I) = (-l)#^ e/ l' 1<s >. 

We denote by T the matrix whose rows and columns are indexed by the a- 
subsets of [—6, a — 1] and the entries of which are: 

T[I J] = l £s{I)Ps if/u W = ( J + 1 ) u {4, (2) 
' 1 otherwise. 

Figure [21 shows two instances of this matrix. The first corresponds to the 
paths with set of steps S — {0, ±1} (Motzkin paths); the second corresponds to 
S = {0, ±1, ±2} (this is known as the basketball problem, and studied in [T]). 
In both cases, we prefer to represent the graph Q, the vertices of which are the 
a-subsets of [—6, a — 1] and the adjacency matrix of which is T. 





Figure 2: The graph Q for the sets of steps S = {0, ±1} (above) and S — 
{0, ±1, ±2} (below). In each case, the subset / of [—6, a — 1] corresponding to 
each vertex is represented by a sequence of • and o (e.g. o»o» corresponds to 
the subset { — 1, 1} of [—2, 1]). The vertex corresponding to the set Io is colored 
gray. 

To help the reader familiarise himself with this definition, we start by stating 
two basic properties of the matrix T. We make use of them later. 

Lemma 2. Let I and J be a-subsets of [—6, a — 1] and assume that T[J, J] ^ 0. 
The two following implications hold. 

1. If —b is in I . then J is such that I U {a} = (J + 1) U {— b}. In this case, 
we have T[7, J] = 

2. If a — 1 is not in J, then I is equal to J + 1. In this case, we have 
T[/,J] = (-l) a /3 a . 
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We can check this lemma by looking at Figure [2J all vertices with a label 
starting with a • have only one outgoing arc, with label /3_ &! all vertices with a 
label ending with a o have only one ingoing arc, with label (— l) a /3 a . 

Proof. By the definition of the matrix T, if / and J are such that T[I, J] ^ 0, 
we have I U {a} = (J + 1) U {s} for some s in [—6, a — 1]. 

Let us prove the first implication. As J is a subset of [—6, a — 1], the point 
—b cannot be in J + 1. Therefore, if — fe is in /, we have s = —6. To prove 
the second implication, we note that the point a is obviously in / U {a} and 
therefore in ( J + 1) U {s}. Therefore, if a — 1 is not in J, we have s = a. □ 

We now define some particular vertices of the graph Q, defined as cyclic 
permutations of Jo- Specifically, if m is such that —b < m < a, let I m be the 
vertex of Q defined as the set [m, m + a — 1] modulo a + b, where the values 
modulo a + b are taken in [—6, a — 1]. Since /_b = I a , these vertices number 
a + b. 

Lemma 3. If < m < a, any walk in the graph Q going from I m visits Iq. If 
—b < m < 0, any walk in the graph Q going backwards from I m visits Iq. 

Proof. The first result stems from the fact that ~b is in I m if < m < a. 
Lemma [2] entails that the only arc of Q going from /,„ goes to I m -i- Thus, any 
walk going from I m reaches Iq in m steps. 

Likewise, the second results stems from the fact that a — 1 is not in I m if 
— b < m < 0. Lemma [5] entails that the only arc of Q going to I m goes from 
I m +i- This means that any walk going backwards from I m reaches Iq. □ 

Proposition 4. The vectors satisfy for k > 0: 

F fe+1 = TF fc . 

Proof. To prove the proposition, we show that for every a-subset / of [—6, a — 1] , 
we have: 

F fc+1 [7] = ^T[7,J]F fe [J]. 

j 

To do this, we let a be a /-permutation of order k + 1 such that [3(a) ^ 0. Let 
s = ct(O) be the head of a and r be its tail, defined above. Since the number of 
inversions created by in the permutation a is equal to the number of elements 
of / lower than s, we have: 

p(a) = s s (I)[3 s (3(t). 

It remains to show that r is a J-permutation, where / U {a} = ( J + 1) U {s}. 
Since f3 s ^ 0, we have s € [—6, a]. Thus, we may write: 

a: N — > {s} U (/ U {a} \ {s}) U N> a+ i. 

From this, we deduce: 

r: N -> JUN>„. 

From its definition, the set J is a priori an a-subset of \—b — 1, a — 1]; however, 
since /3(t) ^ 0, the point —6—1 cannot be in J. This finishes the proof. □ 
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From Proposition 21 we derive the following theorem, which already appears 
in [3] but is found using completely different methods. 

Theorem 5. The generating function of the polynomials Fk is: 

^ k ~ D(z)> 

k>0 y ' 

where D{z) is the determinant of 1 — zT and N(z) is the (Jo, Jo) cof actor of the 
same matrix. 

Moreover, the degree of the polynomial D(z) is ( a+b ) ; the degree of N{z) is 

("•;/•) " 

Proof. The only possible /-permutation of order is the identity, which is an 
Jo-permutation. This implies that the polynomial Fq is 1 if I = Jo and 
otherwise. With Proposition 21 we deduce that the polynomial Ff. is equal to 
the entry (Jo, Jo) in the matrix T fe . 

The generating function J2k>o ^kZ k is therefore equal to the entry (Jo, Jo) 
in the matrix (1 — zT) -1 . The announced form follows from Cramer's rule. 

To compute the degree of the polynomial D{z), we let d — ( a ^ 6 )- Since T 
is a d X d matrix, the polynomial D(z) has degree at most d; the coefficient 
of z d in this polynomial is (up to a sign) det(T). Denoting by &d the set of 
permutations of the set of a-subsets of [—6, a — 1], we have: 

det(T)= Y, eWn^.^]- 

Treed / 

Let f be a permutation with a nonzero contribution in this sum. Lemma [5] 
asserts that: 

1. if -b is in J, then J U {a} = (tt(J) + l) U {-6}; 

2. if a — 1 is not in 7r(J), then J = ir(I) + 1. 

Condition 1 determines tt(I) if —6 is in J; Condition 2 determines tt(J) if —6 
is not in J. The permutation 7r is thus uniquely determined. Moreover, the 
values of T[J,7r(J)] are given by Lemma [21 This gives, up to a sign, the value 
ofdet(T): 

det(T) = ±/3_ 6 ^ a - x } Pa K a \ 

which is nonzero. Therefore, the polynomial D(z) has degree d. 

To compute the degree of N(z), we adopt a more combinatorial point of view: 
we regard the determinant D(z) as the generating function of configurations of 
cycles of the graph Q, counted up to a sign, where z accounts for the number of 
vertices visited by the configuration. Let ttq be the unique permutation, defined 
above, that contributes to the dominant term of D(z); the permutation 7To can 
be interpreted as the only configuration of cycles that visits all vertices of Q. 

The cofactor N(z) is the generating function of the configurations of cycles 
that avoid the vertex Jo. By Lemma [31 such a configuration cannot visit any of 
the vertices I m (since if it would, it would also visit Jo). This means that N(z) 
has degree at most d — a — b. Moreover, we easily check that the vertices J m 
form a cycle of the configuration tto; by removing this cycle, we thus obtain a 
configuration of cycles visiting d — a — b vertices and not visiting Jo. Therefore, 
the degree of N(z) is exactly d — a — b. □ 
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As examples, let us compute the polynomials D(z) and N(z) corresponding 
to the two examples of Figured In the case of Motzkin paths S — {0, ±1}, let 
us take the weights = and lji = w_i = t (i.e., the case of standard Dyck 
paths). Theorem [5] shows that the polynomials Fk satisfy: 

k>0 

which is equivalent the recurrence relation: 

F =F 1 = 1, 

F k = F k -i - t 2 F k _ 2 if k > 2. 

These polynomials are commonly known as the Fibonacci polynomials, due to 
the similarities with the recurrence relation of the Fibonacci numbers. 

In the basketball case S — {0,±1,±2}, let us take the weights ljq = 0, 
u)\ = u)—\ = t\ and u> 2 = = ^i- In this case, the polynomials D(z) and N(z) 
are: 

D(z) = (1 + t 2 zf(l -z- 2t 2 z + tfz 2 + 2t 2 z 2 + 2tjz 2 - tjz 3 - 2t\z z + t\z% 
N(z) = (l + t 2 z)(l-t 2 z). 

As we can see, a simplification by a factor of 1 + t 2 z occurs in the computation 
of the fraction N(z)/ D(z). This implies that the polynomials F k follow a linear 
recurrence relation of order 5 instead of 6, as shown in [T]. The factorisation of 
the polynomial D(z) is also predicted in [3]; it is linked to the fact that the set 
S and the weights /3 S are symmetric. 



3 Bounded meanders 

We enumerate bounded meanders in the same manner as bounded excursions. 
Specifically, the generating function Ek t e counts walks from to £ in the graph 
described by the matrix A k . Therefore, we have 

Ek,e = ~ft^-, (3) 

where F ki e is the (£,0) cofactor of the matrix 1 — Obviously, we have 

Fk,o = Fk. 

We again compute the cofactors Fk t using a transfer matrix method. We 
start by writing Fk/ in terms of the determinant of the matrix 1 — Ak with the 
£th row and the 0th column cut. This reads: 



F M = (-l)'det 



if i < I 
iii>£ 



0<i,j<k-l 
t—1 k-1 

a£G k i=0 i=e 
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For every integer s, we set /3 S = —/3 s +i. We also define, for every permutation 
a of order k, the quantity: 

1-1 k-l 

@l(a) = e(o-) Y[ Pi-<r(i) JJ ft-o-(i)- 

i=0 i=« 

We rewrite the above formula as: 

Let now I be an o - 1-subset of the set [—6 — 1, a — 2]. Let ©1 be the set 
of /-permutations of order k with respect to the set S = S — 1. We define the 
following polynomial: 

We also_denote by Fk,e the vector whose /-component is F^J. 

Let T be the transfer matrix defined in Section [3] corresponding to the set 
S and the weights /3 S . Let U be the matrix whose rows are indexed by the 
a — 1-subsets of [—6 — 1, a — 2], whose columns are indexed by the a-subsets of 
J— 6, a — 1] and whose entries are: 

v[i,J] = \ 1 ^{-^ (4) 

I otherwise. 

We also let Q be the graph with adjacency matrix T. Let TL be the graph 
Q U G, with additional arcs between the vertices of Q and Q coded by the matrix 
U. Examples are shown in Figure [3J 

Proposition 6. The vectors F^^ satisfy for k > £ > 0: 

Ffc,o = UFfe. 

Proof. The proof of the first identity follows the same lines as that of Proposi- 
tion!?] Let J be an a - 1-subset of \—b — 1, a — 2] and a be in ©jS^. Let s be 
the head and r be the tail of a. We have: 

ft+i(<r) = e s {I%[3i{r). 

The rest of the proof is identical to that of Proposition 2J 

To prove jthe second identity, we let I be an a — 1-subset of \—b — 1, a — 2] 
and a be in ©1^. By definition, we have: 

Po(a) = (3(a). 

Moreover, assume that (3q(o-) = /3(c) ^ 0. This means that —6 — 1 cannot be 
in /. Therefore, the set J = I U {a — 1} is an a-subset of [— b, a — 1]. Again by 
definition, the mapping a is in ©j. . This completes the proof. □ 
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Theorem 7. The bivariate generating function of the polynomials F kt i satisfies: 



V F ke u e z k = TF k (u)z k = J* M , 
k>e>o k>o D(uz)D(z) 

where D(z) is the determinant of 1 — zT , D(z) is the determinant of 1 — zT, 
and N(u, z) may be computed as: 

N(u,z)= Cof[I,J ](l -uzT) Cof[io, J](l-zT). 

U[/,J]=1 

The polynomial D(z) has degree { a ^ b ), the polynomial D(z) has degree ("^J), 
and the polynomial N(u, z) has a dominant term in 

( Q+b+1 ) -a- 6-1 ( a+ ?)-a 

Proof. As the only Jo-permutation of order is the identity, Proposition^ entails 
that the polynomial F^t is equal to the entry (Jo, Jo) m the matrix t^UT*^. 
In other terms, Fk t is the generating function of walks from Jo to Iq in the 
graph TL taking t steps in the graph Q and k — £ steps in the graph Q. The 
generating function ^ fe t Fk^u z k is thus equal to the entry (Iq, Jo) in the matrix 
(1 — uizT) _1 U(l — zT) -1 . This yields the announced form. 

Let us now compute the degrees; let d = . Theorem [S] shows that the 
polynomial D(z) has degree d. 

To compute the dominant term of the polynomial N(u,z), we first remark 
that if J and J are such that U[J, J] = 1, we have — b- 1 £ I and a — 1 e J. This 
implies that Lemma[5]is still valid in the graph %. Let J m , for —6—1 < m < a— 1, 
be the vertices of Q defined in the same way as the vertices I m . Since Lemma [2] 
is valid, Lemma [3] is also valid regarding both the vertices J m and I m . 

We now consider the polynomial N(u, z). This polynomial is the generating 
function of configurations in the graph T-L composed of elementary cycles and 
a self-avoiding walk going from Jo to Jo; the variable u takes into account the 
number of arcs of Q in the configuration. By Lemma [31 such a configuration 
cannot visit a vertex J m with < m < a — 1 (since it would contain an arc going 
into Jo), nor can it visit a vertex I m with —6 < m < (since it would contain 
an arc going from Jo). This proves that the dominant term is at most in 



Let 7Tp and ttq be the only configurations of cycles visiting all vertices of the 
graphs Q and Q, respectively. Consider the configuration consisting of: 

• all cycles of 7?o except the one containing the vertices J m ; 

• all cycles of ttq except the one containing the vertices J m ; 

• the self-avoiding walk Jo — > ■ ■ ■ — > I-t — > I a — > ■ ■ • — > lo- 

We check that this configuration contains d — a arcs of Q and d — b — 1 arcs of 
Q . We thus derive the dominant term of the polynomial N(u, z). □ 
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We now consider the generating function Mk of all meanders of height at 
most k regardless of final height. From ©, we find 



M k = (5) 

where the polynomial Gk is the sum of Fk,t for all I between and k. The 
generating function of the polynomials Gk is found by setting w = 1 in the 
expression of Theorem [7J this proves that the polynomials Gk follow a linear 
recurrence relation of order ( a + b ) + = ( a+b+1 ). 

Let us now take the two examples detailed in Section [2] The case of Dyck 
paths (Figure [3J left) is very simple since the graph Q has only one vertex. If 
we set uq — and uj\ = = t, the generating function of the polynomials 
Fk,i is: 

F k ,iu l z k = {l _ tuz){l _ z + t 2 z 2y 
k>e>o y A ; 

In other words, we have: 

Fk,t — t l Fk-t- 

Let us now examine the case where S = {0, ±1, ±2}, ujq = 0, u>i = w-i = t\ 
and 0J2 = u)-2 = t 2 (Figure[31 right). The polynomials D{z) and N(u,z) are: 

D(z) = l-hz- t 2 z 2 - htjz 3 + 4z 4 ; 
N(u, z) = (1 + t 2 z)(l - t 2 z + ti< 2 w^ 2 - 4u 2 z 3 + t\u 2 z A ). 

Once again, a simplification by a factor of 1 + t 2 z occurs in the computation of 
the generating function ^ fc e Fk i eu i z k . 



4 Symmetric set of steps 

We now consider the special case where the set of steps S is symmetric, that 
is, where —S — S and = cjj for all i 6 5 1 . Bousquet-Melou already consid- 
ers this case and shows [3] that the generating function E(t) is canceled by a 
polynomial of degree 2 a instead of ( 2 ^) . 

While we were not able to recover Bousquet-Melou's results with our meth- 
ods, we show that another phenomenon occurs: namely, the polynomials Fk 
factor into two parts, and simplifications occur in the computation of the gen- 
erating functions of meanders Ek,e and Mk- These simplifications are basically 
due to the fact that as S is symmetric, the entry of the matrix Ak is 

identical to the entry (fc — i, k — j). Define the following two matrices: 

\ k J 1 ) / 0<i,i<fe/2 

I UJj-i — Ulk-j-i 
\ /0<i,j<fc/2 

(all values of to s are taken to be if s is not in S). If k is an odd number, both 
matrices have the same dimension and the condition j = k/2 never occurs; if k 
is an even number, the dimension of A^ +1 is one more than that of A^ +1 . In 
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Figure 3: The graphs Q for the sets of steps S — {0, ±1} (above) and S — 
{0, ±1,±2} (below). Graphical conventions are identical to those of Figure 
with the vertex corresponding to Iq colored gray. Arcs with weight 1, coded by 
the matrix U and leading to vertices of the graph Q shown in Figure [5] are also 
shown. 




Figure 4: The graphs corresponding to the three matrices A3, A% and A 3 with 
S = {±1} and uj\ = w_i — t. They differ only by the vertex 3. 
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both cases, the sum of the two dimensions is k + 1. The graphs with adjacency 
matrices A~£ and are illustrated in Figure SJ 

We denote by and FZ the determinants of the matrices 1 — A^_ l and 
1 — Afc_ t , respectively. We also denote by F^ g and F k ~ t the (£, 0) cofactors of 
the matrices 1 — ^4^ and 1 — A^ , respectively. 

Theorem 8. The polynomial F k satisfies for all integers k: 

F k = F k F k- 

Moreover, the generating functions Ek,i satisfy, for all £ such that < £ < k/2: 



E, 



Ek,k-e 



F k,e 



F, 



E, 



k.e 



Ek,k-e 



k,e 



fc+i 



F 



k+l 



Finally, if k is even and £ — k/2, we have 

Ek,e 



F + 



F + 

r k+l 

Proof. Let B = (b 0l ...,bk) be the canonical basis of the underlying vector 
space of the matrix A^- We denote by Bq, B\ and £>2 the following collections 
of vectors: 

Bq = (bi)o<i < k/2] 

&i = (bk/2) if k is even and otherwise; 

£>2 = (bk-i)o<i<k/2- 

Since the entry (i, j) of the matrix A/, is equal to the entry (k — i,k — j), the 
matrix Ak written as a block matrix with respect to the basis (Bo, £>i, S2) looks 
like: 





( * 


u 




A k = 


V 


X 


V 




\ c 


U 


B ) 



Let P be the following passage matrix: 





( 1 





1 \ 


p = 





1 







V 1 





- 1 / 



where 1 denotes the identity matrix of the appropriate dimension. We compute: 

\ 





( B + C 


u 





P- 1 A k P = 


2V 


X 







K 





B-C 













All identities of the theorem are readily derived from this form. 



□ 
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This result shows that the generating functions Ek,t can be computed using 
smaller polynomials than expected, since the matrices and (and thus 
their determinants and cofactors) are about twice as small as A k . Moreover, 
we can derive from the proposition a simplified expression for the generating 
function of meanders Mfc. By writing M% as the sum of E k: i for all < I < k 
and grouping the terms by pairs, we find: 

This expression involves smaller polynomials that ([5]) (a simplification occurs 
by a factor of F k+ 1 ). 

Theorem 9. The generating functions of the polynomials F k and F k are of 
the form: 

^ k ~ D(z 2 ) ' ^ k ~ D(z 2 ) ' 

k>0 y ' k>0 y ' 

where D(z) is the polynomial defined in Theorem^ and N + (z) and N~(z) are 
polynomials. Moreover, the bivariate generating functions of the polynomials 
F^ and F^ e are of the form: 

N+(u,z) , ^ jb N-(u,z) 



>>!>o D(uz*)D(z*) t ^„ ■ D(u^)D(z*) 

where D{z) is the polynomial defined in Theorem^ and N + (u, z) and N~(u, z) 
are polynomials. 

Proof. This theorem is a consequence of a fact hinted at in FigureUJ the matrices 
A~£ and A^ are nearly identical to the matrices A k + and A k - , respectively, 
where k + = If J and k" = I ^j^J ■ More precisely, the term uik-j-i is zero if 
k — j — i > a; since j < fc/2, this is true whenever i < fc/2 — a. 

For simplicity, we only prove the results associated to the matrix A k , but 
the case of A k is identical. The proof follows the same techniques used in the 
previous sections. If / is a a-subset of [—a, a — 1], we denote by F^" +1 the vector 
whose /-coefficient is: 

Assume now that fc — 1 > 2a and examine the first term of the product: the 
above remark entails that Wfc_ CT (i)_i = 0, which means that the first term is 
equal to (3 a (i)-i- This allows us to repeat the proof of Proposition 0) As the 
matrix A k minus its first row and first column is equal to A k _ 2 , we find: 

In the same way, we define the vectors e ; by repeating the proof of Pro- 
posotionEl we find, if fc is sufficiently large: 

F M = UF,. 

All the identities of the theorem are derived from these recurrence relations. □ 
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As an example, we take the Fibonacci polynomials, corresponding to the set 
S = {±1} (see Section [2]). The values of the polynomials F^ and FZ are given 
by: 



F+=F k -tF k ^; 
F 2k = F k + tFk-i; 



F + 



2fe+l — 



Fk+i — t 2 Fk-i; 



2fe+l — 



F k . 



Obviously, these four sequences of polynomials follow the same recurrence rela- 
tion as the polynomials F k . 
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